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Previously suggested correlations of BL Lac type objects with the arrival directions of the ultra-high- 
energy cosmic ray primaries are tested by making use of the HiRes stereoscopic data. The results of the study 
support the conclusion that BL Lacs may be the cosmic ray sources and suggest the presence of a small (a 
few percent) fraction of neutral primaries at _B > 10^^ eV. 
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1. MOTIVATION 

With the exception of a small number of highest- 
energy events, the bulk of ultra-high energy cosmic rays 
(UHECR) is most naturally explained by acceleration in 
active galaxies [2 121 El- This explanation is supported 
by clustering of UHECR observed in the AGASA and 
■Yakutsk datasets |31E1I1] ^^ and by correlations between 
arrival directions of UHECR and a particular class of ac- 
tive galactic nuclei, the BL Lacertae objects (BL Lacs) 
[niiniini (see also ReL HU). 

The correlations with BL Lacs were observed in 
AGASA and Yakutsk datasets. To test them, inde- 
pendent datasets are needed. Such a dataset has been 
■recently published by the HiRes collaboration 'S'. It 
'consists of 271 events with reconstructed energy E > 
10^^ eV which were observed in the stereoscopic mode. 
The energies of the events are not published; this makes 
reconstruction 12 of arrival directions in the Galactic 
magnetic field impossible. Despite the small exposure of 
HiRes in the stereo mode as compared to AGASA, the 
excellent angular resolution of 0.6° makes HiRes data 
quite competitive (or even superior) in testing correla- 
tions with BL Lacs under the assumption of neutral pri- 
mary particles. This test is the purpose of the present 
work. More precisely, we take three previously studied 
sets of BL Lacs without any further selection and test 
them for correlations with the entire set of 271 HiRes 
events. 



2. FLUCTUATIONS OF CORRELATION 
SIGNAL 

At the current low statistics the correlation signal is 
not expected to be robust "7. To illustrate this impor- 
tant point consider the following example. Assume that 
a given set of sources in a given exposure time produces 
5 events in average. Let the average number of rays 
coincident with sources by chance be 2. The average 
number of correlating rays observed in such an experi- 
ment would therefore be 7. Both the real events from 
sources and background chance coincidences are Pois- 
son distributed. Thus, in a particular realization the 
number of correlating rays would vary roughly within 
7 ± v^- It is therefore likely that two identical exper- 
iments will see 10 and 4 correlating rays, respectively. 
The first one would see a strong correlation signal (10 
events at 2 expected with the hypothesis of no correla- 
tions; probability p ~ 5 x 10~^), the second one would 
see a signal compatible with the uniform background (4 
events at 2 expected, p ~ 14%). 

It is worth stressing that, in this situation, the final 
conclusion would be that the correlation exists, because 
the first of the two experiments is incompatible with 
uniform background (while both are compatible with 
the signal). Thus, at present small statistics one should 
expect that some of the experiments see no correlations. 
Such fluctuations indeed exist, e.g., in the autocorrela- 
tion signal ^ 8^ . 



^) e-mail: gorby@ms2.inr.ac.ru 

^' HiRes collaboration does not observe significant cluster- 
ing 13 |S| but, at the present level of statistics, this does not 
contradict to the previous results IHIllUI . 



3. BL LAC SAMPLES 



There were three samples of BL Lacs discussed in the 
literature, all three drown from the catalog of quasars 
and active galactic nuclei |15| : 
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(1) The set of 22 most powerful BL Lacs, Ref. [TT] . 
This set was found by adjustment of cuts on magnitude 
and radio flux in the catalog of conflrmed BL Lacs using 
the requirement of the best correlation signal with the 
most clustered subset of AGASA and Yakutsk cosmic 
ray events. Neutral primaries were assumed. 

(2) The set of 14 BL Lacs which contains poten- 
tially 7-ray-loud objects, Ref. jEl- This subset was ob- 
tained by the positional cross-correlation of confirmed 
BL Lacs with the catalog of EGRET sources ITBl. Both 
assumptions of neutral and charged primary particles 
were tested using the most clustered subset of AGASA 
and Yakutsk events. 

(3) The set of all confirmed BL Lacs with the single 
cut on visual magnitude, mag < 18, Ref. 12 . This 
set of BL Lacs contains 156 objects. Both neutral and 
charged primaries were tested using the whole available 
set of AGASA events, E > 4 x 10^^ eV. Better corre- 
lations were found under the assumption of positively 
charged primary particles. 

Note that the subsets (1-3) are not independent. BL 
Lacs which correlate with the UHECR in the set (1) 
are also present in the set (2). The cut on magnitude, 
mag < 18, being the single cut for the set (3), is also 
imposed in the set (1) among other cuts. 



4. PROCEDURE 

Our analysis is based on the calculation of the an- 
gular correlation function by means of the algorithm 
described in Ref. 11 . The statistical significance of 
correlation is estimated by testing the hypothesis that 
the UHECR and BL Lacs are uncorrelated. This is done 
as follows. For a given set of sources and the angle 6, 
we count the number of pairs "cosmic ray - source" 
separated by the angular distance less or equal S, thus 
obtaining the "data count" . We then replace the real 
data by a randomly generated Monte-Carlo set of cos- 
mic rays and calculate the number of pairs in the same 
way, thus obtaining the "Monte-Carlo count" . We re- 
peat the latter procedure many times noting cases when 
the Monte-Carlo count equals or exceeds the data count. 
The number of such cases divided by the total number 
of tries gives the probability P{S) which characterizes 
the significance of correlations at a given angular scale 
S. The smaller is this probability, the stronger (more 
significant) is the correlation. 

The Monte-Carlo events are drown from the 
isotropic distribution according to the acceptance of a 
given UHECR detector. For the HiRes experiment in 
the stereo mode we use the zenith-angle and azimuth- 



angle distributions given in Ref. il7i| and the sidereal 
time distribution published in Ref. jl8| . 

The probability P{5) depends on the choice of S. 
The optimal value of S is clearly the one where the ex- 
pected signal is strongest (e.g., expected Pth(^) is low- 
est). In the case of point sources and neutral primary 
particles the shape of Pth{S) can be calculated by the 
following Monte-Carlo simulation. We start by gener- 
ating a cosmic ray set which is artificially correlated 
with BL Lacs. This is achieved by replacing a given 
(small) number of events in a random set by the events 
correlated with randomly chosen BL Lacs. To generate 
an event correlated with the BL Lac we assume the 2- 
dimensional Gaussian distribution normalized in such a 
way that the circle of the radius 0.6° around the BL Lac 
contains 68% of events. The resulting correlated set is 
then treated as the real data, i.e., the probability P{5) 
is calculated. The bin-by-bin median values of the prob- 
abilities obtained for many different sets give Pth{S)- In 
the real data, following Ref. ^\, we measure correla- 
tions at that value of d where Pth(<5) is minimal. 



5. RESULTS 

In Figures ^ - 13 we present the expected (thin line) 
and observed (thick line) P{S) for the HiRes dataset of 
UHECR and three samples of BL Lac objects described 
in Section U. There are no significant correlations with 
the subsets of 22 and 14 BL Lacs. In both cases only 
two BL Lacs contribute to data counts within 1 degree. 
In the subset of 22 BL Lacs those are TXS 0751-^485 
and OT 465 (the latter object was among 5 BL Lacs 
which contribute to correlations in Ref. jE]). In the 
subset of 14 BL Lacs two contributing objects are PKS 
1604-^159 and lES 1959-^650. It is worth mentioning 
that lES 1959-1-650 is a confirmed TeV source and con- 
tributes to correlations in Ref |13j . 

The set of 156 BL Lacs with mag < 18 shows 
rather strong correlation. We list in Table 1 the cos- 
mic rays and BL Lac's separated by less than 1°. The 
expected probability Pth(<5) has minimum at S — 0.8. 
At this value of S the real data give the probability 
P(0.8) = 4 X lO"**. The data count is 11, while the 
Monte-Carlo expectation due to random uncorrelated 
background is ~ 3. This corresponds to about 8 events 
from sources in the overall HiRes data-sample of 271 
events. Thus, a ^ 3% fraction of neutral particles in 
the total cosmic ray flux may be sufficient to explain this 
correlation. This is consistent with a fraction of clus- 
tered events observed by AGASA at these energies [201 ■ 
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FIG. 1. P{5) for the set of 22 BL Lacs and HiRes 
stereo events. The thick line shows data, the thin line 
shows Pth{S) obtained in the Monte-Carlo simulations 
in which 2 events are from the BL Lac sources. 




FIG. 3. P{5) for the set of 156 BL Lacs and HiRes 
stereo events. The thick line shows data, the thin line 
shows Pth(<5) obtained in the Monte-Gar lo simulations 
in which 9 events are from the BL Lac sources. 




FIG. 2. P{5) for the set of 14 BL Lacs and HiRes 
stereo events. The thick line shows data, the thin line 
shows Pth{S) obtained in the Monte-Garlo simulations 
in which 2 events are from the BL Lac sources. 



6. CONCLUSIONS 



Cosmic ray 


BL Lac 


a, deg. 


S, deg. 


name 


z 


17.8 


-12.5 


RBS 0161 


0.234 


48.5 


5.8 


RX J03143-f 0620 


? 


118.7 


48.1 


TXS 0751-^485 


7 


123.8 


57.0 


RX J08163-f 5739 


? 


137.2 


33.5 


Ton 1015 


0.354 


162.6 


49.2 


MS 10507+4946 


0.140 


169.3 


25.9 


RX J11176-h2548 


0.360 


209.9 


59.7 


RX J13598+5911 


7 


226.5 


56.5 


SBS 1508+561 


7 


229.0 


56.4 


SBS 1508+561 


7 


253.7 


39.8 


RGB J1652+403 


7 


265.3 


46.7 


OT465 


7 


300.2 


65.1 


lES 1959+650 


0.047 



Table 1. Cosmic rays and BL Lac's from the sample of 
156 objects separated by less than 1°. Note that SBS 
1508+561 correlates with two rays. 



We have assumed — as a "null-hypothesis" — that 
BL Lacs and the new HiRes stereo data are uncorre- 
lated. We have tested this hypothesis for 3 particular 
subsets of BL Lacs studied previously and found an ex- 
cess of pairs "BL Lac - cosmic ray" in one of them, 
which would occur with the probability 4 x 10""' for a 
random set of UHECR. If these sets were independent, 
the final significance would be the best probability found 
multiplied by the number of tries. Since the samples are 



not independent, the penalty factor is smaller. In any 
case, the null-hypothesis is rejected with the probability 
^ 10"'^. Our analysis thus confirms the association of 
UHECR with BL Lacs. 

Only one of the three examined subsets of BL Lacs 
— the largest one — shows significant correlations. This 
may indicate that those BL Lacs which are sources, have 
similar and small luminosity in UHECR. This conclu- 
sion is also consistent with the statistics of the clus- 
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tering of UHECR events, which suggests few hundred 
sources currently contributing to the UHECR flux, see 

Refs.EHEHUni- 

The observation of correlations at angles which are 
much smaller than the typical deflection of a charged 
UHECR particle in the Galactic magnetic field suggests 
that (at least) a few percent of UHECR are neutral. 
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